Spectral Function in Mott Insulating Surfaces 
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We show theoretically the fingerprints of short-range spiral magnetic correlations in the photoe- 
mission spectra of the Mott insulating ground states realized in the y/3 x \/3 triangular silicon 
^v^j . surfaces K/Si(lll)-B and SiC(OOOl). The calculated spectra present low energy features of mag- 

netic origin with a reduced dispersion ~ 10 — 40 meV compared with the center-of-mass spectra 
bandwidth ~ 0.2 — 0.3 eV. Remarkably, we find that the quasiparticle signal survives only around 
the magnetic Goldstone modes. Our findings would position these silicon surfaces as new candidates 
to investigate non-conventional quasiparticle excitations. 

It is well known that the strong correlation in electronic systems with an odd number of electrons per unit cell 
may give rise to Mott insulating (MI) ground states, in contradiction with band theory predictions. Generally, the 
presence of metal-transition ions with unpaired electrons are responsible for the MI properties of many compounds 
like the undoped high-T c cuprates. However, it has also been shown that surfaces characterized by a y/3 x y/3 
arrangement of silicon dangling-bond surface orbitals can provide the ideal conditions to realize the MI phenomena 
^ without metal-transition ions[l], 0, 0> l1- Although the on-site Coulomb potential U is not so large, the surface state 
bandwidth is reduced due to reconstructions, becoming the correlation effect important. Recently, Weitering et aZ.pJ 
used moment-resolved direct and inverse photoemission spectroscopy to demonstrate that the triangular interface 
K/Si(lll)-(v3 x Vo)-B (hereafter K/Si-B) has a MI ground state. Even though such an insulating state is not 
i-^j [ necessarily unstable toward antiferromagnetism due to the lack of perfect nesting, it has also been argued that K/Si- 
B is the first experimental realization of a triangular Heisenberg spin-1/2 model with a 120° Neel order. However, 
standard density-functional methods combined with exact diagonalization studies suggest a Mott insulator with a non- 
magnetic ground stateQ. In addition, it was speculated that the reconstructed triangular surface SiC(0001)-(\/3 x \/3) 
(hereafter SiC) also presents a 120° Neel ordered ground stateQ. Since direct measurements of magnetic order in 
£Nj , such surfaces are difficult to implement, the magnetic properties of these MI ground states still remain unanswered^. 
On the other hand, nowadays it is possible to perform higher resolution photoemission experiments than the previous 
ones, which would allow a detailed analysis of the single-hole properties in these MI ground states. These issues 
motivated us to study theoretically the single-hole dynamics in a triangular antiferromagnet (AF) with two objectives: 
i) to investigate the effect of a frustrated magnetic order in the quasiparticle (QP) behavior, and ii) to obtain 
spectroscopic fingerprints of magnetic order through the photoemission spectra calculated for realistic parameters of 
^vq ' the surfaces K/Si-B and SiC. A careful comparison with future higher resolution spectroscopy experiments could give 
, some information about the underlying short-range magnetic structure. 

As a consequence of the magnetic order, we find a strong k-dependence of the spectral function structures. In 
particular, we show the emergence of spectral weight of magnetic origin between the Fermi level and the measured 
surface state band, with a strongly reduced dispersion. For the surfaces K/Si-B and SiC, our results show a remarkable 
vanishing of the QP weight for a large region of the Brillouin zone (BZ) outside a neighborhood of the magnetic 
Goldstone modes. 

To tackle this problem we assume the t—J model, which is supposed to capture the low-energy physics involved 
, in photoemission spectroscopy of antiferromagnetic Mott insulators [fj. In order to solve the model, we use the self- 
consistent Born approximation (SCBA) complemented with exact diagonalization studies. For non-frustrated AF, it 
.£h | has already been established the success of the SCBA to reproduce exact diagonalization results on small clusters Q 
■ and angle-resolved photoemission spectroscopy (ARPES) experiments^. 

We assume a surface ground state with a magnetic wave vector Q = (47r/3,0) lying in the surface plane x— z, and 
spin waves as the low energy excitations. We perform a unitary transformation to local quantization axis so as to 
have a ferromagnetic ground state in the z' direction. Then, we use the spinless fermion = h\, of, = hiS^ , and 

the Holstein-Primakov Sf ~ h{ a i + a i)> &i ~ f ( a l — a i)i "Sf — \ ~ a \ a i representations. These are replaced in 
the t—J model keeping the relevant terms up to third order. After a lengthy but straightforward calculation, the 
Hamiltonian results 
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In the Hamiltonian l|TJl. £k = — ^7k and cj q = | J^/ (1 — 37 q )(l + 67 q ) are the hole and magnon dispersions, re- 
spectively. 7k = J2s cos(k.<5) and /3k = J2s sin(k.<5) are geometric factors, with 5 the positive vectors to nearest 

neighbors, and k varying in the first BZ of the V3 x V3 surface. The bare hole-magnon vertex interaction is 

i 

defined by Mk q = i(/?kV_ q — Ac-qWq) with the Bogoliubov coefficients it q = [(f + 37 q /2 + (x> q )/2w q ] 2 and u q = 

- 

sign(7 q ) [(I + 7 q /2 — 3w q )/2aj q ] 2 . The free hopping hole term implies a finite probability of the hole to move without 
emission or absorption of magnons because of the underlying non-collinear magnetic structure. The hole-magnon 
vertex interaction adds another mechanism for the charge carrier motion which is magnon-assisted. The latter is the 
responsible for the spin-polaron formation when a hole is injected in a non- frustrated antiferromagnet Q ■ We will 
show the existence of a subtle interference between both processes that turns out to be dependent on the momenta. 
An important quantity that allows us to study the interplay between such processes is the retarded hole Green func- 
tion that is defined as G£(oj) = (AF\h-k ^ +i ^+_ H ^ h^ k \AF), where \AF) is the undoped antiferromagnetic ground state 
with a 120° Neel order. In the SCBA the self energy at zero temperature results 
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We have solved numerically this self-consistent equation for S^fu), and calculated the hole spectra function A^(u>) = 
— ilmGkfy) and the quasiparticle (QP) weight Zk = f 1 — ^"%7^) U k > where the QP energy is given by E k = 

Sk(i?k). Before we discuss the results it is important to mention that, unlike previous works 0, we will concentrate 
on the behavior of the photoemission spectra for realistic parameters of the surfaces K/Si-B and SiC, what implies a 
strong coupling regime and a careful extrapolation to the thermodynamic limit (we have studied cluster sizes up to 
2700 sites). 

For the K/Si-B (SiC) surface, photoemission spectra give a bandwidth W ~ 0.3 eV (~ 0.2 eV) for the occupied 
surface electronic band0, 0- These experiments together with inverse photoemission indicate an effective on-site 
Coulomb repulsion U ~ I — 2 eV (~ 1.5 — 2 eV). Electronic band calculations give a similar value for U and a wider 
surface bandwidth W ~ 0.61 eV (~ 0.35 eV), due to the neglection of correlation effects, leading to a positive t ~ 0.07 
eV (~ 0.04 eV) 0,H| . The large U/W ratios indicate the presence of strong electronic correlations and a considerable 
suppression of charge fluctuations. In addition, scanning tunneling microscopy measurements on the SiC surface 
show a clear evidence of a Mott insulating state ^(j- Theoretical works on the Hubbard model suggest that these 
surfaces would be located in the antiferromagnetic region of the phase diagram^] • All these estimations point out 
to J/t ~ 0.1 — 0.4 for both surfaces. In that parameter range, we have observed negligible quantitative changes in our 
results, so we take J/t = 0.4 as a reference value. 




FIG. 1: Spectral functions versus frequency for J/t — 0.4 and N — 21 corresponding to the momenta F = (0,0), A = 
|y( — 1, 3v3) and B = |y (2, \/3) shown as filled circles in the inset of the left panel (the crosses represent the other momenta). 
The solid and dashed lines are the exact and SCBA results, respectively. 



In order to test the validity of the SCBA we have compared its results for the single-hole dynamics with Lanczos 
exact diagonalization (ED) calculations on small clusters. In general, we have found a good agreement for all momenta 
of the BZ and for J/t < 1. In Fig. ^we show the spectral functions for a cluster size of N = 21 and in the strong 
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FIG. 2: Spectral functions for realistic parameters of the silicon surfaces (J/t = 0.4) calculated at three different points of the 
BZ (see inset Fig. [3J. Notice that the Fermi level has been placed on the right. 



coupling regime J/t — 0.4, for three momenta of the BZ. There is a very good agreement between ED and SCBA 
results in the whole range of energies. This concordance gives a strong support to our analytical approach and so, 
from now on, we will focus on the SCBA spectra in the thermodynamic limit. 

In Fig. [3 we show the spectral function structures, which can be traced back to the distinct hole motion processes 
mentioned above. The spectra extend over a frequency range of ~ 9t, that is, the non-interacting electronic bandwidth. 
Similar values has been found in the experiments on the silicon surfaces (HQ. For momentum T (upper panel) both 
processes contribute coherently to the quasiparticle excitation at the top of the spectra. There is also a small incoherent 
component centered around 6t below the QP peak. When moving away from T, there is a spectral weight transfer from 
the low energy coherent sector to higher energies. At M' (middle panel) the coherence is completely lost, surviving 
only a magnetic tail at low energy and a structureless background. At X (lower panel), corresponding to the magnetic 
vector Q, the coherence emerges once again, giving rise to the QP. In the background there is a broad resonance 
related to the free hopping hole motion with a finite lifetime of order ~ 4 J. 

In Fig. |3| we show the intensity of the QP peak along high symmetry axes and, in the inset, the region of the 
BZ where the quasiparticle weight is finite for J/t = 0.4. Besides the strong k-dependence, it can be observed that 
there is no QP for momenta outside the neighborhood of the magnetic Goldstone modes (k = 0, Q). The existence 
of the quasiparticle peak is also observed for an ample range of J/t values. In particular, as J/t increases the shaded 
areas in the inset of Fig. [3J where the quasiparticle weight is finite, increase and only for values J/t > 2.5 there is 
quasiparticle peak all over the Brillouin zone. The non existence of quasiparticle is a striking manifestation of the 
strong interference between the free and magnon-assisted hopping processes. It is expected that this feature could be 
observed in any system with short-range non-collinear correlations. 

The available photoemission spectroscopy data for the silicon surfaces have a low energy resolution (~ 100 — 200 
meV) 0,12 and it is impossible to discern the energy structure of the surface spectra. The surface bands obtained 
experimentally are correctly reproduced by the center-of-mass of our spectra. The latter coincides with the free 
hopping dispersion (see Fig. as it is expected from the exact treatment of the first spectral moment of our theory. 
Using the transfer integral t obtained from first-principles local-density- approximation calculations, Q our bandwidths 
compare very well with the experimental ones, reflecting the narrowing induced by the electronic correlation. For 
the K/Si-B (SiC) surface we have obtained W ~ 0.3 eV (~ 0.18 eV). In Fig. 0]we show the center-of-mass and the 
photothreshold dispersions for the SiC surface. When the quasiparticle exists, the photothreshold corresponds to the 
QP energy excitation. As it can be observed, the underlying magnetic structure changes the center-of-mass minimum 
at the K point to a QP energy local maximum, nearly degenerate with the hole ground state momentum Y. There is 
also an appreciable reduced bandwidth ~ 10 meV of the photothreshold dispersion in comparison with the measured 
surface state bandwidth ~ 200 meV. Whereas for the K/Si-B surface the values are ~ 40 meV and ~ 300 meV, 
respectively. 



4 




FIG. 3: QP intensity along the T — K — M — T path for J/t = 0.4. Inset: 
QP weight is finite. 



\/3 x \/3 Brillouin zone. In the shaded areas the 




FIG. 4: Surface band structure for realistic parameters ( J/t=0.1, t = 0.04 eV) of the SiC(0001)-(\/3 x \/3) surface. The dotted 
line is the photothreshold energy, the solid one is the QP dispersion, and the dashed one is the center-of-mass spectra band. 
Energies are given relative to the measured bulk valence-band maximum (VBM). Experimentally the Fermi level is located at 
2.3 eV above the VBM 



In conclusion, we have studied theoretically the hole spectral function in the triangular t— J model for realistic 
parameters relevant for the silicon surfaces SiC(0001)-(\/3 x \/3) and K/Si(lll)-\/3 x \/3- Assuming the presence 
of a long-range magnetic Neel order we have observed the emergence of low energy features of magnetic origin with 
a reduced dispersion band. As the photocmission spectrum is not sensitive to the asymptotic low energy magnetic 
properties of the system, we speculate however that it could give important information about the presence of short- 
range magnetic order. We have also obtained an unexpected vanishing of the QP weight for a large region of the 
Brillouin zone for these Mott insulating surfaces. Our theoretical predictions could provide useful ground to the 
analysis of future improved photocmission experiments. Using a simple and reliable analytical method (SCBA), we 
have found clear signatures of interesting physics caused by strong electronic correlation on simple silicon surfaces. 

This work was done under PICT Grant No. N03-03833 (ANPCyT) and was partially supported by Fundacion 
Antorchas. 
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